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A quantum gravity-gradiometer consists of two spatially 
separated ensembles of atoms interrogated by pulses of a com- 
mon laser beam. Laser pulses cause the probability ampli- 
tudes of atomic ground-state hyperfine levels to interfere, pro- 
ducing two motion-sensitive phase shifts which allow the mea- 
surement of the average acceleration of each ensemble, and, 
via simple differencing, of the acceleration gradient. Here I 
propose entangling the quantum states of the two ensembles 
prior to the pulse sequence, and show that entanglement en- 
codes their relative acceleration in a single interference phase 
which can be measured directly, with no need for differencing. 

PACS numbers: 39.20.-Fq, 03.75.Dg, 03.67.-a, 03.67.Lx 



1. Introduction 

Inertial sensors based on matter-wave interferometry 
are a rapidly developing technology [1 — 4]. One class of 
such sensors are gradiometers designed to measure linear 
acceleration gradients, typically gradients caused by in- 
homogeneities in the gravitational potential. The state- 
of-the-art design Q for an atom-wave gravity gradiome- 
ter consists of a pair of atom- wave accelerometers sepa- 
rated spatially by a fixed distance and direction. Each 
of the two accelerometers contains an ensemble of laser- 
cooled atoms, through which a common pair of (counter- 
propagating) Raman laser beams are pulsed in a carefully 
controlled sequence to drive Rabi oscillations between 
atomic ground-state hyperfine levels. The sequence and 
timing of the laser pulses are adjusted so that the hyper- 
fine ground and excited states of any single atom interfere 
during the atom's motion through the laser field. At the 
end of the sequence, the probability of finding an atom in 
its excited hyperfine state is a simple function of a rela- 
tive phase, which accumulates between successive pulses 
in proportion to the atom's acceleration (relative to the 
inertial frame defined by the laser beam) averaged over 
its fiight time. Subsequently, the resulting pair of ob- 
servations of this phase (one for each ensemble) can be 
subtracted to obtain a measurement of the relative ac- 
celeration, and, upon dividing the result by the separa- 
tion length, a measurement of the acceleration gradient. 
I will give a more quantitative review of this quantum 
interferometry technique in Sect. 2 below. A detailed de- 
scription of the experimental setup can be found in ^] 
(see esp. Fig. 1 there). 

Even though subtracting the two acceleration measure- 
ments allows a number of correlated noise sources to can- 



cel out as common modes y| , it is still desirable to avoid 
such differencing of two nearly-equal measured quanti- 
ties which are corrupted by (uncorrelated) noise. In this 
letter, I will show that by prior-entangling (prior to the 
laser pulsing sequence) the quantum states of atoms in 
the two ensembles (such that every entangled pair has 
one atom in each ensemble) it becomes possible to en- 
code the relative-acceleration information directly in a 
single interference phase, thereby eliminating the need 
for the differencing of two separate phase measurements. 
This idea is inspired by the recent discovery of a sim- 
ilar quantum algorithm which uses entangled states to 
synchronize atomic clocks non-locally . The Quantum 
Clock Synchronization (QCS) algorithm relies on prepar- 
ing the atoms "stabilizing" each atomic clock in a special 
entangled quantum state whose time evolution reduces to 
a pure multiplicative phase as long as each atom in the 
pair evolves under the same unitary transformation. By 
contrast, the Entangled Quantum Interferometry (EQI) 
algorithm to be described below relies on transforming 
the pair's quantum state into a sensor which senses the 
difference between unitary evolutions of the entangled 
atoms. Put another way, the EQI algorithm is a spa- 
tial analogue of the QCS algorithm, in exactly the same 
sense as atomic clocks are temporal analogues of quan- 
tum (atom- wave) interferometers Q. 

2. Overview of matter- v^rave gravity gradiometry 

Consider an atom in its internal ground state |0), mov- 
ing through the inertial frame defined by the Raman 
beams of an atom interferometer. [Here "inertial frame" 
has the same meaning as in general relativity: in the pres- 
ence of a gravitational field, my coordinate system will 
always be locally-Minkowski (free- falling, or geodesic- 
normal coordinates) so that the gravitational potential 
and its first derivatives vanish at the origin.] The wave 
function can be written as a tensor product of the atom's 
internal and external Hilbert-space states: 

-^0:p{x,t) = |0) (g) j a(fc,u;)e*(^-^-'^*) Sk 

= |0)®i^(f,t) , (1) 

where a{k,uj) defines the wave packet representing the 
external (space-time) part ipi^^t) of the quantum state, 
uj is related to k via the usual dispersion relation (plus 
an additive constant corresponding to the atom's internal 
energy level) , and p denotes the average (in general time- 
dependent) momentum: 



p= -ifiVlV') = (27^)^ y hk\a{k,uj)\'^ d?k . (2) 

The need for this cumbersome subscript notation 5'o:p 
will become clear in a moment. Since the two Raman 
beams are detuned in frequency by an amount which 
corresponds to the energy difference between the ground 
and excited states |0) and they might stimulate the 
ground-state atom to absorb a photon from one beam 
and emit a lower frequency photon into the other, result- 
ing in a transition to the excited state This stimu- 
lated transition gives the atom a net momentum kick in 
the amount h{ki — ^2) = hK^ where fci, k2 {\ki\ > \k2\) 
denote the wave vectors corresponding to the (counter- 
propagating) Raman beams. [As the vectors ki, i = 1, 2, 
point in opposite directions, the magnitude of K is typi- 
cally twice that of the fci . Note that fl is the total energy 
absorbed by the atom, which includes not only the tran- 
sition energy (within the appropriate bandwidth), but 
also the recoil contribution.] The wave function of the 
same atom excited in this way from its ground state can 
then be decomposed, similarly to Eq. (1), in the form 

"^i-.n+ni} i^^t) = Il> ® / a(fc,c.)e^[(^^+^)-^~-("+")*l d'k 



(3) 



= \l)®^Pix,t) e'(K-S-m 



Since the wave packet 'ip{x,t) is common to both the 
ground and excited state wave functions, I will suppress 
it (as an overall "normalization" factor) in what follows. 
Note that whenever an atom interrogated in the interfer- 
ometer is excited from its ground state |0 : p) to the inter- 
nal excited state |1) by Raman pulses, it always picks up 
an associated excess energy-momentum (K, fl) as speci- 
fied in Eq. (3); in other words, the excited atom's total 
quantum state becomes 1 1 : p + KK) . Conversely, a stim- 
ulated transition in the reverse direction results in the 
transformation \1 : p) — > \0 : p — KK). This correla- 
tion between an atom's internal state and its momentum 
during the laser interrogation process is the key feature 
enabling atom interferometry |Q . 

In the interferometer, atoms are manipulated by laser 
pulses of two basic kinds: A Hadamard pulse iJ| (a 7r/2 
pulse followed by the spin operator cr^), whose action is 



H^\0:p) = -j={\Q:p) + \l:p + hK)) , 
H^\l:p + hK) = ^{\Q:p)-\l: p + HK)) , (4) 

and a "double" Hadamard pulse (a tt pulse followed 
by the spin operator cr^), whose action is 



H,\0:p) 
H^\l:p + hK) 



\l:p + hK) 
|0:p) . 



Consider a single atom's flight across the laser held dur- 
ing the interrogation phase: Initially, at times t < ti, say, 
the atom is in the ground state |0 : p). At time t = ti the 
first HiL pulse hits, and the atom's state becomes (up to 
an overall phase factor which I will always ignore) : 



vl'(ti) = i= (|0:p) + |l:p + fti?)) 



(6) 



Since atoms are well-localized spatially, their wave func- 
tions "^{x, t) are sharply peaked around their average po- 
sition (5'(i)| X l^'(i)) at all times t. When an atom is ex- 
cited by Raman pulses to a superposition of its internal 
states, its distribution becomes bimodal, but still sharply 
peaked around the two positions (modes). Therefore, to 
a very good approximation, I can follow the evolution of 
the atom's wave function during its flight in the inter- 
ferometer by examining it in the vicinity of the atom's 
(average) position(s) as a function of time t. (See |^ for 
a more rigorous treatment.) Accordingly, at any time t 
after ti and before the next Raman pulse hits, the state 
Eq. (6) evolves, according to Eqs. (l)-(3), as 



1 
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\0:p)+ e^-^-(^'-^i)-''"(*-^*i)|l : p + KK) 

ti<t<t2, 



(7) 



where xt = x{t), xi = x{ti) etc. is shorthand notation for 
the atom's position at the different times, and an overall 
phase factor as well as the common wave packet ip{x,t) 
are suppressed as advertised. At a time t ~ t2 > ti, a,n 
pulse hits, and the atom's new state becomes 



^{t2 
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\l:p + hK) 



_|_ gi-ff ■(?2-Sl)-ii^(t2-tl) |Q . 



(8) 



At any time t after t2 and before the next pulse, the state 
Eq. (8) evolves as 



1 



^iK-(xt—X2)-in{t — t2) 



1 : p + hK) 



gl/?-(x2-Xl)-iO(t2-tl)|Q . p^ 

Vt : t2<t<t3 . 



(9) 



(5) 



At time t ^ t^, the second and last iJ| pulse hits, and 
transforms ^'(/;) into the state 

^{t3) = i \^e''^■(^=>-^2)-^n(ts-t2)QQ : p) _ |i : p-^ f^x)) 

^^K■{x2-x,)-^n{t2-t^)^^Q : p) + |1 : p+ hK))] . (10) 



So as a result of the pulse sequence — iJ^ — , an 
atom which enters the interferometer in the ground state 



|0 : p) ends up, at the end of the last Hjl pulse (at t = tj,)^ 
in the state Eq. (10) which can be rewritten as 



*(i3 



|0:p) 



+ (e«=^i _e*e32) \l:p + hK)) 



where for i,j — 1, 2, 3, . . . 

9y = K ■ {xi - Xj) - fl{ti - tj) . 



(11) 



(12) 



After absorbing a factor of e'®^^ into the overall phase 
multiplying '^(ta), Eq. (11) takes the more familiar form 



*(i3 



(1 



\0:p) 



(1-e'®) \l:p + hK)) 



(13) 



where Q = 832 — 021- The phase 8 can now be observed 
by measuring the relative abundance of ground vs. ex- 
cited state atoms in the state ^'(^3). Specifically, the 
fraction Pi of excited-state atoms in the state Eq. (13) is 
given by 



e 

sm- I - 



(14) 



and, similarly, the fraction Pq of ground state atoms is 
Po = |1 e^Q|2/4 = cos2(e/2). On the other hand, 
substituting Eqs. (12) in 6 = 632 — 621 gives 

Q = K- [x{t3)~2x{t2) + x{ti)] ~ n{t3-2t2+ti) . (15) 

If the Raman pulses making up the sequence Hzl — — 
H2L are aligned in time such that t2 = ti + T and ^3 = 
t2+T for some common interrogation time T, the second 
term in parenthesis in Eq. (15) vanishes, and, using the 
standard Taylor expansion of x{t), we can rewrite the 
first term in the form 



K ■ x{t2) + 0{T^ 



(16) 



Here the magnitude of the 0{T^) remainder is down by 
a factor of order T/Ta relative to the first (acceleration) 
term, where Ta ^ ||a||/||a|| denotes the timescale over 
which the atom's acceleration a{t) = x(t) varies. There- 
fore, as long as the time scale Ta is much larger (as is 
typically the case) than the interrogation time T, mea- 
surement of the phase O following the pulse sequence 
-H^— iJ^ — iJ| yields a direct measurement of the 
inertial acceleration component K ■ a. 
3. Interferometry with entangled atoms 

For a gradient measurement, ordinarily it is necessary 
to apply two separate but simultaneous measurements 
of acceleration to two different ensembles of atoms sepa- 
rated spatially by a fixed distance. By using atoms whose 



internal states are pairwise entangled in a distributed 
fashion, however, it is possible to obtain the acceleration 
gradient information with just one phase measurement 
as I will show now: 

Start, at time t ~ ti, with atoms in the two ensembles 
A and B in the pairwise entangled initial state 



*(ii 



35(1° 



pa) A |1 : PS + hK)B 



V2 

|1 ■.pA + hK)A \Q.Pb)b 



(17) 



In other words, in some subset (which my measurements 
will be directed at) of ensemble A, every atom is en- 
tangled with some atom in ensemble _B, and, at some 
initial time t — ti, each entangled pair is (upto an over- 
all phase) in the state Eq. (17). Focus attention now on 
a single such entangled pair. Since the state Eq. (17) is 
an energy eigenstate of the total internal Hamiltonian, 
its explicit time evolution consists of a pure multiplica- 
tive phase only. Implicitly, it evolves in time kinemati- 
cally with the inertial motion of the atoms; according to 
Eqs. (l)-(3), this evolution takes the form 



*(0 



V2 



e'^-^"^^^-^^'^ |0 : pa) A |1 : Pb + hK)B 
\1:PA + hK)A |0 : pb)b 

yt: h<t<t2, 



(18) 



where XAt = XAit), xbi = 2^5(^1), etc., and an over- 
all phase factor as well as the common wave packet 
ipAixA,t)'4'B{xB,t) are suppressed as in Eq. (7) above. 
For increased readability, let me introduce the following 
abbreviation which I will use throughout the rest of this 
letter: for E = A, B, 



\0)e = \0:pe)i 



\1)e = \l:pE + hK)i 



(19) 



Now suppose that, at a time t = t2 > ti, a, common Hj^ 
pulse is applied to both atoms A and B (simultaneously) 
in the state Eq. (18). Then the new state of the pair at 
t = t2 becomes 



*(t2 



1 
72 



^iK-{xB2-XBl} 



1)a |0)i 



^iK-{xA2-SAl) 



\0)a |1) 



(20) 



At any time t after t2 and before the next pulse, the state 
Eq. (20) evolves as 



JK-{xAt-^A2+XB2-3Bl) 



|1)a |0)s 

iK-(xBt-XB2+3A2-XAl) |Q^^ \\) B 



yt: to <t<U 



(21) 



At a time t — > t2, & final, common Hjl pulse is applied 
to the entangled pair A and B, and transforms the state 
Eq. (21) into 



1 



2^/2 



,-t(*-*32 + *^2l) 



(|0)a-|i)a)(|0)b + |i>b) 



— e 



)a + |i>a)(|0)j 



, (22) 



where, for i, j = 1, 2, 3, 



K-{xAt-XAj), ^^ij=K ■ {xbi - xbj) ■ (23) 



After absorbing a factor of e*^* into the overall 

phase multiplying 5* (is), Eq. (22) takes the more man- 
ageable form 



1 



e^*(|OU-|l)^)(|0)B + |l)B) 



2^2 

(|0U + |i)a)(|0)b-|i)b) 



(24) 



with the phase observable $ given by 

$ = $^32 + $^21 - $'^32 - ^''^21 

= K ■ (Afa - 2 Af2 + Axi) , 



(25) 



where Ax = — x_b and Axi = x^ii — XBi for «,j 
1,2,3,... Collecting terms in Eq. (24), 



"fits) 



1 



(e**-l)|0)A|0)B + (e** + 1)10)^11)5 

(e'*-l)|l)^|l)B] . (26) 



2\/2 

(e^* + l)|l)^|0) 



The phase $ can now be observed by measuring the rela- 
tive abundance in the state \I'(t3) of those pairs in which 
both atoms are in their excited states. [Even though 
^(^3) is not stationary, its time evolution consists of pure 
phases multiplying each of the coefficients in Eq. (26), so 
the observation time is not critical.] Specifically, accord- 
ing to Eq. (26), the fraction Pn of those pairs in which 
both the j4-atom and the B-atom are in the excited state 
|1) is given by 



Pi 



11 



1 



1 



1 



- sm — 
2 \ 2 



(27) 



On the other hand, if the Raman pulses making up the 
sequence —H,^ — Hjl are aligned in time as before so that 
t2 — ti+T and ^3 = t2+T for some common interrogation 
time T, we can write, using Eq. (25), 



'S^ ^ K ■ Ax{t2) + 0{T^) 



(28) 



where the magnitude of the O(T^) remainder is again 
down by a factor of order T/Ta relative to the first (ac- 
celeration) term, with Ta denoting the timescale over 



which the atoms' acceleration varies. As long as the time 
scale Ta is much larger than the interrogation time T, 
measurement of the single phase observable $ following 
the pulse sequence — i?^ — 7J| on the entangled pair 
Eq. (17) yields a direct measurement of the inertial ac- 
celeration gradient K ■ Aa = K ■ {Ha — as) between the 
ensembles A and B. Signal sensitivity of $ to the accel- 
eration gradient Aa can be enhanced using tricks similar 
to the one discussed in the note [|| . 
4. Producing the entangled atom pairs 

A number of methods for entangling two-level atoms 
have been discussed in the literature and some imple- 
mented in experiments. For completeness, I will sketch 
two such methods which can be used to produce states 
of the form Eq. (17). 

One method, originally proposed in Cabrillo et. al. , 
involves pumping two spatially separated atoms A and 
B into their excited states, so that the joint system is 
initially in the product state 

\^)a\1)b. 

A single-photon detector, which cannot (even in princi- 
ple) distinguish the atom from which a detected photon 
arrives, is placed halfway between the atoms. Such a 
detection scheme can be implemented either by using a 
detector which is inherently insensitive to the direction of 
the photon that triggers it, or by using two detectors sit- 
uated at the output ports of a symmetric beam-splitter 
on which the light paths from the two atoms are sym- 
metrically impingent, thus erasing the which-path infor- 
mation. When one of the atoms makes a transition to its 
ground state, and the detector registers the emitted pho- 
ton, the result of its measurement is to put the combined 
two-atom system into the entangled state 

^ (|0>^|l)B + e^1l)A|0)B) , (29) 



where the phase is a measure of the relative delay be- 
tween the paths reaching the detectors from A versus 
from B, and can be adjusted to have any desired value 
[e.g. as in Eq. (17)] by adjusting the position of the de- 
tector relative to the atoms. 

A second method, recently investigated in detail by 
Haroche and coworkers |l^ is (in very rough outline) the 
following: Start with a single-mode cavity whose exci- 
tation frequency is tuned to 17. Send the pair of atoms 
A and B into the cavity one after the other, with atom 
B first. Initially, both atoms and the cavity are in their 
ground states: 



|o)a |o)s |o) 



EM 



(30) 



where |0)£;a/ denotes the vacuum state of the cavity. Af- 
ter atom B is in the cavity, apply a 7r/2-pulse on it, which 
transforms the state Eq. (30) into 



A 



— \Q)^{\0)s\1)em-\1)b\0)em) 



(31) 



When both atoms are in the cavity, apply a second, tt- 
pulse, this time on the atom A, thereby transforming the 
state Eq. (31) into 



1 



{\1)a\0)b - |0)a|1)s)®|0) 



EM 



(32) 



which, for the atom-pair A and B, is in the desired 
form Eq. (17) up to an overaU phase factor. Consider- 
ing that a number of proof-of-principle experiments con- 
firming it have already been performed successfully , 
this method (or, more precisely, a suitable variant which 
maintains the necessary correlation between atomic en- 
ergy levels and linear momentum as discussed Sect. 2 
above) is likely to be the method of choice in practi- 
cal implementations of the entangled-atom gradiometry 
algorithm. 

5. Quantifying the improvement 

To what extent does the entangled-interferometry tech- 
nique proposed in Sect. 3 represent an improvement over 
the standard, "classical" atom-wave gradiometry as re- 
viewed in Sect. 2? To answer this question quantitatively, 
I will now carry out a simple analysis of the improve- 
ment in the signal-to-noise ratio expected with entangled- 
interferometry. 

I will take the differential phase-shift as my observed 
signal, so in a gradient measurement with both classical 
and entangled interferometry, the signal is the dimen- 
sionless quantity 



S = K-AdT^ , 



(33) 



where a = x, and other symbols are defined as in 
Eqs. (16) and (28), with the same approximation that 
allows neglecting terms of order 0{T^) being assumed 
implicitly. In this discussion, for ease of notation I will 
label the two atomic ensembles hy i — 1, 2 as opposed to 
by A and B before (in Sects. 2 — 3). In a "classical" gra- 
dient measurement, two ensembles containing M atoms 
each are used to obtain the two acceleration measure- 
ments, which are later differenced to obtain the signal 



^ = ei - 62 , 



(34) 



where the phase observables Qi are related to the ob- 
served abundance (fiuorescence) of excited-state atoms 
in each ensemble i = 1, 2 via Eq. (14): 



M 



(35) 



Here Ali denotes the number of excited-state atoms at the 
end of the interrogation process in ensemble i, i — 1, 2. 



For comparison, imagine a similar experiment using 
the entangled algorithm of Sect. 3, involving M pairs 
of entangled atoms interrogated in the pulse sequence 
of Eqs. (17) — (23) with the same interrogation-time con- 
stant T. At the end of that sequence, the number of atom 
pairs Mil in which both atoms end up in their excited 
states is measured, which yields the signal 



5" = $ 

where according to Eq. (27) 
1 



■ sm 



Mil 
M 



(36) 



(37) 



For simplicity, I will now make the additional assump- 
tion that the only sources of noise in both kinds of sig- 
nal measurements (classical and entangled) are (i) phase 
noise, and (ii) shot noise due to the statistical fluctua- 
tions in the observed numbers of excited atoms: in Mi 
and M2 for the classical case, and in Afn for the entan- 
gled case. I will assume a binomial distribution for these 
random variables with M trials and with probabilities 
Pi, P2, and Pii, respectively. The expected RMS shot 
noise in each number measurement Mi is then given by 



AM, = y <M,2> - <Af,>2 = ^P,{1~P,)M , (38) 
and the expected RMS noise in Mn is, similarly. 



AA/n = \J <Mii^> - <Afn>2 = ^Pii{l-Pii)M . 

(39) 

In the classical measurement [Eqs. (34) — (35)], the shot 
noise contribution to each Qi measurement can be writ- 
ten as 



dM, 



AM,, 



1, 2 



(40) 



Since the two measurements of Qi are completely inde- 
pendent, the total variance (noise squared) in the quan- 
tity S — AQ is the sum of the variances in each O^; 



Nc = {Ni^ + N2') ~' « V2N1 = V2 -T^AMi , (41) 

UlVli 



where the approximate equality holds because of the as- 
sumption that the signal <C 1. Substituting Eqs. (35), 
(38) and (14) in Eq. (41) gives 



(42) 



In the entangled measurement, according to Eqs. (36) — 
(37) 



5" = $ = 2 arcsm ^ / sa 2v2 



M 

and the shot noise contribution is 
'2(1 -Pii) 



oMii 



M 



(43) 



where the second equahty foUows from Eq. (39), and the 
third by the assumption that the signal <C 1. Re- 
markably, the shot noise contributions in the classical 
and entangled measurements are identical. Since a typ- 
ical classical gradient measurement operates in a regime 
where the number of atoms M in each ensemble i is large 
enough to make shot noise negligible jj] , the same operat- 
ing conditions will make it negligible also in the entangled 
case, and therefore phase noise is the dominant source of 
error in both kinds of measurement. 

To compare phase-noise performance of the two kinds 
of measurement, I will model the noise as simply propor- 
tional to the magnitude of the fluorescence signal M: 



> - <ej>2= sin^ 



for the classical case, and 

<$2> _ <$>2= 0-21 sin^ 



<e,> 



<<!>> 



1, 2 



(45) 



a gradiometer in low earth orbit (common-mode acceler- 
ations due to atmospheric drag being typically of order 
10~^g or higher) the phase-noise suppression Ne/Nc via 
entangled interferometry can be expected to be a reduced 
but still significant two orders of magnitude or better. 

6. Future work 

The Entangled Quantum Interferometry algorithm as 
presented here can be generalized to design sensors of 
other "distributed" quantities (such as higher derivatives 
of the gravitational curvature tensor); these generaliza- 
tions will be discussed in a future publication. 
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i = 1, 2 (46) 



in the entangled case, where cr is a constant. By the same 
arguments as above [Eqs. (40) — (41)], the phase noise in 
the differential acceleration signal is, for the classical case 



/- . <9i> 
iVc « V 2 (7 sm — - — 



and for the entangled case 
1 



T2 



(47) 



(48) 



Therefore the ratio of phase noise in the entangled mea- 
surement to that in the classical measurement can be 
written in the form [cf. Eqs. (36), (34), and (33)] 



1 sin 



i sm 



<ei> 



1 <$> 

2 <ei> 



I Ac 



(49) 



For a typical measurement performed on earth with a Im 
separation between the two atomic ensembles, the earth's 
~ 3000E gradient field (1E= 10~^sec~^) and Ig gravity 
suggests suppression of phase noise via entanglement by 
a factor of up to 1.5 x 10^^, or more than six orders of 
magnitude. For measurements performed in spacecraft 
in low earth orbit, the ambient acceleration ||a|| will be 
due to non-isolated common-mode vibrations originating 
from drag forces acting on the spacecraft. Therefore, for 
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